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Abstract 

The critical exponents vl2jVL2 o,nd 7^2 of a uniaxial Lifshitz point are cal- 
culated at two-loop level using renormalization group and ei-expansion tech- 
niques. We introduced a new constraint involving the loop momenta along the 
competition axis, which allows to solve the two-loop integrals. The exponent 
7L2 obtained using our method is in good agreement with numerical estimates 
based on Monte Carlo simulations. 
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I. INTRODUCTION 



The Lifshitz point occurs in a variety of physical systems and has been extensively studied 
over the last twenty-five years JIJH. It appears in High-T^ superconductivity polymer 
physics ferroelectric liquid crystals @, |Kj, etc.. In magnetic systems, the uniaxial 

Lifshitz critical behavior can be described by an axially next-nearest-neighbor Ising model 
(ANNNI). It consists of a spin-| system on a cubic lattice (d = 3) with nearest-neighbor 
ferromagnetic couplings and next-nearest-neighbor antiferromagnetic interactions along a 



single lattice axis [ill. The competition gives origin to a modulated phase, in addition to 
the ferromagnetic and paramagnetic ones. In spite of having several modulated phases, it 
was shown recently that around the Lifshitz critical region, a simple field-theoretic setting 



can be defined for this ANNNI model 0]. In general, the antiferromagnetic couplings can 
show up in m directions. In that case, the system possesses the m-fold Lifshitz critical 
point. Here we are going to focus our attention in the uniaxial case (m = 1), since some 
materials present this type of critical behavior. MnP was studied both theoretically and 
experimentally and displays the uniaxial behavior fl3^|Tij[ . 



Theoretical studies involving the uniaxial Lifshitz point have been put forth using analyt- 
ical and numerical tools. Examples of the latter are high-temperature series expansion [16| 
and Monte Carlo simulations |LT|,|rT|| . Conformal invariance calculations in d = 2 (in the 
context of strongly anisotropic criticality) and e-expansion techniques |], [19], |2(| have 
been the main analytical tools available to dealing with this kind of system. 

From the field-theoretic point of view, the critical dimension of a scalar field theory 
describing the uniaxial Lifshitz critical behavior is found (by using the Ginzburg criterion) 
to be d c = 4.5. The expansion parameter is €l = 4.5 — d, where d is the space dimension 
of the system under consideration. (In the pure Ising model the expansion parameter is 
e = 4 — d). As is well known, the critical exponents for the Lifshitz point at one- loop 
approximation have the same dependence in as those from the pure Ising model have in 
e. Of particular importance is the effect of the mixing of the two momenta scales, i. e., along 
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and perpendicular to the competing axis. One can choose a convenient symmetry point to 
fix the external momenta scale in the quartic and quadratic directions. Then, the one-loop 
integral contributing to the four-point function, needed to find out the fixed point, can be 
performed without any approximation. The choice which simplifies the referred integral is to 
set the external momenta scale in the quartic direction to zero. The solution of this integral 
yields a leading singularity and a regular term in e^. The leading singularity in can be 
chosen to be the same as that obtained in e when solving the analogous one-loop integral for 
the Ising model (by absorbing a convenient angular factor in a redefinition of the coupling 
constant), even though the coefficients of the regular terms in and e are slightly different 



and depend on m in the Lifshitz case ||20|| . Thus, although the multiplicative factors to be 
absorbed in the coupling constant are different in the two cases, the critical exponents have 
the same dependence on the expansion parameter. 

We can then ask ourselves if this happens to be true when we proceed to calculate the 
critical exponents in higher orders in the perturbative expansion. Going one step further 
to evaluate higher order integrals in the loop expansion would be highly desirable using 
the same line of reasoning. In this way, all the dependence in the external momenta is 
along the (d — 1) directions, perpendicular to the competition axis. In the resulting e L - 
expansion, the leading singularities can be chosen equal to those coming from the theory 
without competition (see below). The nontrivial new features of this expansion around the 
usual quadratic field theory are the coefficients of the subleading singularities and of the 
regular terms, which are no longer rational numbers. This approach would allow to treat 
this system properly in a perturbative expansion at least for correlation functions along the 
perpendicular directions to the competition axis. A better comprehension of this procedure 
might shed light in the perturbative study of higher order derivative field theories. In this 
sense, the Lifshitz critical behavior seems to be the natural laboratory to study higher order 
field theories in a perturbative framework. 

We report on what we believe to be the first study of critical exponents at two-loop 
order for the uniaxial Lifshitz point. Using A0 4 field theory and the expansion in powers of 



eL = 4.5 — d in the Lifshitz critical point, we give a detailed account of the calculation of the 
exponents z/^ 2 and r]L2, which by now can be viewed as a worked out example of a recent 



generalization obtained for anisotropic behaviors ||20|| . In order to solve higher- loop integrals 
we introduce a constraint relating the loop momenta in internal and external subdiagrams 
along the competing axis. The results for these integrals are consistent with homogeneity of 
the Feynman integrals in the external quadratic momenta scale. The exponents Vl2 and r)L2 
are associated with the directions perpendicular to the competition axis. (The exponents 
vii and T]L4: associated with the competition axis are not going to be considered here and we 
shall analyse them in another work.) We then obtain the exponent 7^2 via scaling relations. 
The paper is organized as follows. In section 2 we introduce the notation and calculate 
the relevant integrals for the determination of the critical exponents at two-loop level. We 
present the critical exponents rj L2 , vl2 and 7^2 in section 3. In section 4 we discuss our results 
and compare the 7^2 exponent with numerical estimates based on Monte Carlo methods and 
high-temperature series. 



II. CALCULATION OF HIGHER LOOP INTEGRALS 

The Lifshitz critical behavior can be described using a modified A0 4 field theory. The 
bare Lagrangian associated with the uniaxial critical behavior is given by: 

L = \\ V? 0o T + \\ V(d-i) 0o T + 5o^| Vi 0o T + ^o0o + ^o0o- (!) 

The competition is responsible for the appearance of the quartic term in the free prop- 
agator. The Lifshitz critical point is characterized by the values to = 5q = 0. From now 
on, this is the case which interests us in this work. First, we are going to compute the 
renormalized coupling constant at the fixed point. In order to find out the factor that shall 
be absorbed in the coupling constant, we quickly review the one- loop contribution to the 
four point function fl2| . The relevant integral is: 

d d ~ l qdk 
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(2) 



The external momenta are k! along the quartic (competing) direction and p along the 
(d — 1) quadratic directions. We then choose a symmetry point which simplifies the integral 
at external momenta k! = 0, p 2 = 1. Using Schwinger's parameterization we get 

f d d ~^qdk f°° f°° ( f°° \ 

J WTb+mw+¥) = J° 'o daida xL * ex p(-(«i+^)j 

x J d d ~ x q exp(-(«i + a 2 )q 2 - 2a 2 q.p - a 2 p 2 ). (3) 
The q integral is straightforward. It can be performed to give 



J d d 1 q exp(— (oti + a 2 )q 2 - 2a 2 q.p - a 2 p 2 



= 7^-ir(— — )(oti + a 2 ) 2 exp( — ). (4) 

z z a.\ + a 2 

The next step is to compute the k integration, which is: 

roo | 1 | 

2 dk exp(-(ai + a 2 )k A ) = -(<*i + a 2 )~T(-). (5) 
jo z 4 

Replacing equations © into equation @ together with the value p 2 = 1, one finds 

d d ~ l qdk \ 1 1 

( fc 4 + (g + p )2 )(fc 4 + g2)y | p2=i 4 I 2 J l 4 J 



OO /"OO 



x / / da\da 2 exp( — )(ai + a 2 ) ("V+D. (6) 

JO jo «i + «2 

We can perform one of the integrals in the Schwinger parameters using a change of 
variables. Set v = " 2 , and u — a\V. The integral over it can be done, and we are left 
with 



/ / (foida^ exp( ){oti + a 2 ) K 2 

jo jo ai + a 2 

= r(2 - + -)) Z 1 d«(«(i - ^))^)- 2 . (7) 



2 4 y jo 

Now we make the continuation d = 4.5 — e^. One obtains a result in terms of Gamma 
functions with non integer arguments. A useful identity involving the expansion of Gamma 
functions around a small number is given by: 

T(a + bx) = T(a) [ 1 + bx ip{a) + 0{x 2 ) 1 , (8) 



where ip(z) = 4-lr\F(z). This allows one to obtain the e^-expansion when the Gamma 
functions have non integer arguments. Replacing Eqs. (7), (8) into Eq. (6), we obtain: 

I 2 = ir(7/4)r(l/4)^_ 1 l(l + i 2 e L ), (9) 

where i 2 = 1 + \(ip(l) — i/j(1))- We absorb in the coupling constant a geometric angular 
factor, which is |r(7/4)r(l/4)5 , d _i, where S d = [tf-^T^)]- 1 0. Then the redefined 
integral is: 

(10) 



ir(7/4)r(i/4)5 d _i 



or 



~h = -{l + i2e L ). (11) 

We suppress the tilde hereafter to simplify the notation. We have to keep in mind that 
we should divide out this factor for each loop integration. We now turn our attention to 
higher-loop integrals. In practice, we have to calculate the two-loop integrals I^sp = h, 
^h\sp = I3 and the three-loop integral ^h\sp = ^5, m order to find the fixed point at 
two-loop level and the critical exponents. The subscript SP is used to denote our choice of 
the subtraction point. They are given by (see Figure 1): 



d d 1 q 1 d d 1 q 2 dk 1 dk 2 ^ 



{ql + kf) (q 2 2 + k$) (( qi + q 2 + P y + (k x + k 2 + k'Y) 
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d d - l q x d d - Y q 2 d d - Y q z dk x dk 2 dk z 



(ql + kf) (ql + ki) (ql + **) (( Ql +q 2 - pf + (h + k 2 - k'Y) 



(qi + q 3 - PY + (h + h~ k'Y 



(13) 



d 1 qid d v q 2 dkxdk 2 



(qf + kf) ((P - qiY + (K' - hY) (ql + K) 



x 



(qi - q 2 + P3) 2 + (h - k 2 + k' 3 Y 



(14) 



In the first two integrals, p is the external momentum (associated with the two-point 
vertex) along (d — 1) directions, whereas k' is the external momentum along the competi- 
tion axis. Inside the integral I4, P = p\ + P2, with pi,P2,P3 being the external momenta 
(associated with the four-point vertex) along the quadratic directions, and K' = k[ + k' 2 , 
with k' x ,k' 2 ^k' z the external momenta along the quartic direction. The symmetry point is 
chosen as follows. We set all the external momenta at the competition axis equal to zero. 
For the four-point vertex, the external momenta along the quadratic directions are chosen 
as Pi-Pj = j(45jj — 1). We fix the momentum scale of the two-point function through 
p 2 = k 2 = I. 

Now we can study the solution of the higher-loop integrals shown above. Consider the 
integral I 3 . With our choice for the quartic external momenta it is given by: 



r d d 1 q 1 dk 1 r d d 1 q 2 dk 2 

3 ~ J (q 2 i + kj) J (q 2 2 + k\) {{ft, + q 2 + p) 2 + (h + k 2 f) ' 1 J 

First, we perform the integral over the internal bubble, i.e., we integrate over q 2 and ki. 
In order to solve the internal bubble we demand that the loop momenta k\ should be related 
to &2- Note that we could have chosen the other way around, since the integral is symmetric 
under the exchange k\ < — ► k 2 , q\ < — > ?2- There are two issues which need to be emphasized 
here. First, the Lifshitz point condition eliminates the quadratic part of the momenta along 
the competition axis, for 5 = 0. Second, the remaining quartic part of the integral mixes 
the two loop momenta (ki,k 2 ) in different subdiagrams, yielding crossed terms which are 
difficult to integrate. Indeed, using Schwinger parameters and carrying out the integration 
over q 2 first, we obtain 

,d-l s f d d ~ l qidh 



/s( P ,o) = ^ d _ m r(^i)/. 



q\ + (k 2 ) 2 

I" I" da 1 da 2 {a 1 + a 2 )^ 11 exp{ °^-{ qi + p) 2 ) [ dk 2 e- ai{k ^ V a2((fcl+/C2)2)5 \ (16) 

JO JO Ctl + «2 J 



In order to integrater over k2, we have to expand the argument of the last exponential. 
This results in a complicated function of a±, «2, k\ and k 2 , which has no elementary primitive. 



Considering the remaining terms as a damping factor to the integrand, the maximum of the 

integrand will be either at k\ — or at ki = —2k 2 . (The most general choice k x = —ak 2 

yields a hypergeo metric function.) The constraint eliminates the crossed terms and is the 

simplest way to disentangle the two quartic integrals in the loop momenta. The choice 

k\ = —2k 2 implies that k\ varies in the internal bubble, but not in an arbitrary manner. Its 

variation is dominated by k 2 through this constraint, which eliminates the dependence on 

k\ in the internal subdiagram. Integrating over k 2 yields a simple factor to the remaining 

parametric integrals (over the two Schwinger parameters) proportional to (oti + a 2 )~* ■ After 

solving the parametric integrals we find: 

= f d^q.dh I 

3 J (9? + tf)[(?i+p) 2 ]* 2 ' 
We use Schwinger's parameterization again to solve this integral along the quartic direction. 

We obtain: 

Js = 12 I 7^77 — ~ ^2l^ ' ( 18 ) 
J [giH(9i +P) 2 } 2 

The difference with respect to the pure 4 theory is that after performing the quartic integral, 
we get a exponent for the quadratic part of the momenta which is not an integer. At this 
point, one can use Feynman parameters to solve the momentum integrals. The dependence 
of the integral in the external momenta is proportional to (p 2 ) 1 ~ €i , in conformity with the 
homogeneity of the Feynman integrals in the external momenta scale. Deriving with respect 
to p 2 and setting p 2 — 1, we find : 
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3 " 7e L 
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l+[i2 + ^)e L \. (19) 



The integral I4 can be calculated in a similar fashion. First, we set the external momenta 
along the competing direction equal to zero. Therefore 

d d ~ 1 q±dki 



= ( l . 

J (<?? + ki) ((P - qi) 2 + kf) 

d d ~ l q 2 dk 2 



(ql + kl) [(qi-q 2 +pz) 2 + {k 1 + k 2 y 



(20) 



where we changed variables from k 2 — > —k 2 . We set k\ = —2k 2 in the internal bubble q 2 , k 2 
(as we did for 7 3 ), and integrate over q 2 , k 2 . We then have 



d^qxdh 1 

Pz) 2 ]^ 



(21) 



(qf + kf) ((P - ft) 2 + kf) [( gi 
We use Schwinger's parameterization to get rid of the k\ integral. After performing the 
change of variables used to calculate the one-loop integral and a rescaling, we can solve one 
of the parametric integrals to get 

d d ~ 1 q 1 



h = h / dz 

J o J {q 2 - 2z P. qi + zP 2 Y [{q x + p 3 ) 2 } 



(22) 



In order to perform the integral over qi, we make use of a Feynman parameter obtaining 
/, = -/ ! I. - — c(-) J dyy* (l-y)a« 



1 1 -i/ 2 fl-^^)) r( ^ 

2 V 2 V 4V r (^\ Jo 



ke L -l 



x dz 
Jo 



yz(l - yz)P 2 + - y)p\ - 2yz(l - y)p 3 -P 



(23) 



There is a subtlety that needs to be analyzed with care. Here we proceed in complete 
analogy to the pure 4 theory The integral over y is singular at y — 1 when = 0. 
We add and subtract the value of the integrand in the last integral at the point y = 1 



yz(l - yz)P 2 + y(l - y)p\ - 2yz{\ - y)p 3 -P 



-CL 



[z(l-z)P 



2-i-tL 



__ l ]u , [yz{l-yz)P 2 + y(l- y)p 2 3 - 2yz{\- y)p 3 .P] \ | ^ 



(24) 



z{l - z)P 2 

As y — > 1 the logarithm is zero when = 0, leading to a well defined result for the integral 
over y. The coefficient of the integral is proportional to — , which cancels the in front of 
the logarithm. The logarithm contributes only to the order e° L and can be neglected. We 
then find 



2ei 



1 + 3 i 2 e L 



(25) 



Finally, let us describe the computation of the three- loop integral J5. At zero external 
momenta along the competition axis, this integral reads: 
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d d 1 q 1 dk 1 r 



d d - 1 q 2 dk 2 



(q! + kf) J (g 2 2 + k\ 



)((qi + q2-p) 2 + (k 1 + k 2 y) 



x 



/ 



d d 1 q 3 dk 3 



(26) 



(qi + H) (fai + qs - p) 2 + (h + W 



The integral is symmetric in q 2 < > q 3 , k 2 < > k 3 . As the loop momenta are dummy 

variables, the two internal bubbles are really the same. That is why we do not need more 
than one relation among the loop momenta in the external and internal bubbles, even though 
this is a three- loop diagram (with two internal bubbles). As a matter of fact, we can solve 
the two integrals independently in the following way. In the internal bubble ^3, k 3 we set 
k\ = —2k 3 , as well as kx = —2k 2 over the other internal bubble q 2 , k 2 . Apparently we 
have two different relations among the loop momenta, but one of them is artificial. This 
means that the two internal bubbles give the same contribution, i.e. the integration over 
the internal bubbles is proportional to If. Thus 



We integrate first over ki and proceed analogously as before, to find the following result: 



We stress once again that the constraint preserves the physical principle of homogeneity 
in all the higher-loop Feynman integrals in the external momenta scale, which is consistent 
with scaling theory. With the integrals calculated in this way, we can find out the exponents 
as is going to be shown in the next section. 




(27) 




(28) 



III. CRITICAL EXPONENTS 



To compute the critical exponents associated to the ferromagnetic planes at the Lifshitz 
critical point, one may use the standard field-theoretic approach |^TJ. This is possible since 
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no new renormalization constants need to be introduced in this caseQ. From the results 
for 7g, Eq. (|T9|), and Eq. (p8|), we note that these integrals do not have the same 
leading singularities as in pure 4 theory. As an approximation, we introduce a weight 
factor for the two point vertex function, in order to identify the leading singularities with 
the ones appearing in the pure 4 field theory. This factor is | for the integrals above. 
This approximation is suitable when one considers the generalization for the m-fold case 
with m ^ 8. In this way, we have a smooth transition from the Isinglike case (m = 0) to 



the general Lifshitz anisotropic critical behavior (m 7^ 8) p0 |. The bare and renormalized 
quantities are related through O = %1 <fi an d u = Z^ 2 Z u u, where <po and Ao = k el uq are 
the bare parameters in Eq. (1). As usual, = 1 + 5^ and Z u = 1 + 5 U are the wave- function 
and coupling constant renormalization constants, respectively. In addition, we introduce the 
composite field renormalization constant Z^p. , and also Z^p. = Z^p Z^ = 1 + 8^2 . We have 

/91nw N 

7* = A.^ (29) 
9 In 
ou 

The fixed point u* is given by the solution of the equation (3 u {u*) = 0, and the critical 
exponents by the relations r/ L2 = 7</>(w*) and = 2 — r] L2 — 7|(w*)- In case the order 
parameter has a O(N) symmetry, the formulas relating the integrals computed in section 2 
and the renormalization constants defined above are: 



5^ = B 2 u 2 I' 3 + (2B 3 1 2 I' 3 - B 3 I' 5 ) u 3 + 0(u 4 ) 



5 U = 3A 1 ul 2 + 3{&A\ if - A ( 2 ] if - 2A { 2 2) J 4 ) u 2 + 0{u 3 ) (30) 



Ci u I 2 + (C 2 i 2 2 - d i 4 ) u 2 + 0{u 3 



1 This is not valid in the calculation of the critical exponents along the competition axis. 
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where A x = (iV + 8)/18, = {N 2 + 6N + 20)/108, A { 2 2) = (5N + 22)/54, B 2 = (JV + 2)/18, 
5 3 = (JV + 2)(iV + 8)/108, Ci = (N + 2)/6, and C 2 = (N + 2)(N + 8)/36. 

With this information we compute the fixed point at two-loop level. We expand the 
dimensionless bare coupling constant uq in terms of the renormalized coupling u and the 
€l parameter. Using Eqs. ( p9|) and ( |3~0|) , we find the fixed point for the O(N) symmetric 
theory in the following form: 



6 



u 



8 + N 



/ 4(5iV + 22) 
(8 + iV) 2 



(2 + N) 
(8 + N) 2 



Therefore, the exponents r?£ 2 and z/^ 2 are given by: 



VL2 



2 + N 
8 + Nf 



(31) 



(32) 



, ,3 (2 + N) 



+ iV) J 



/ 4(5iV + 22) 
i (8 + iV)2 



1\ 1 _ (2 + TV) 
2 7 2 7 (8 + iV) 2 



(33) 



^L2 



1 1 2 + N 



4 8 + N 



(34) 



+ 



1 (2 + iV) 



2(14iV + 40) z 2 - 2(2 + N) + (8 + JV)(3 + N) 



(35) 



,■ v 3 + Nf 

Now using Fisher's law along directions perpendicular to the competing axis, namely 
7L2 = ^2(2 — T]L2>i the exponent 7^ 2 can be written as: 



1 2 + N 
7L2 = 1 + 2 tL 8+~N 



(36) 



+ 1.(2 + iV) 



12 + 8A" + N 2 + Ai 2 (20 + 7iV) 



(37) 



4 (8 + AO 3 

Previous results in the literature only yielded the exponent 77^2 at 0(e 2 L ) and the exponent 
VL2 at 0(ei)- For this uniaxial case, our results express the critical exponents in a higher 
order in ei compared to earlier investigations. A detailed comparison with other methods 
is provided in the next section. 
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IV. DISCUSSION 



First of all, our result for the exponent 7/1,2 is in agreement with Mukamel's calculation |2! 
at 0(t 2 L ). Therefore, our method is equivalent to integrating over the momentum shell as 
was done in his work using the Landau- Ginzburg- Wilson Hamiltonian approach. 



For the ANNNI model, 7^2 = 1-4 ± 0.06 is the former Monte Carlo output [ IT | , whereas 



the best estimates from the high-temperature series is 7^2 = 1-62 ± 0.12 [IB|. Note that we 
use the subscript 7^2 instead of 75, since it was shown recently that the exponents parallel 
and perpendicular to the competition axis obey independent scaling laws [ 2"3]] . Our two- 
loop calculation obtained from the e^-expansion via the scaling law (when neglecting 0(e\) 
terms) in three dimensions yields 7L2 = 1-45. This agrees (within the error bar) with the 
former Monte Carlo result, the difference being of order 10~ 2 . 

Nevertheless, the most recent high-precision numerical Monte Carlo estimate for the 



ANNNI model yielded 7^2 = 1.36 ± 0.03 JTTj. In order to figure out how to extract the 
best numerical results from the e^-expansion when the parameter is not small (which is 
the case for d — 3), a comparison with the Ising model is worthwhile. For the exponent 
7 in three dimensions, the e-expansion gives a contribution of 0.167 at 0(e) and 0.077 at 
0(e 2 ) |21[. The O(e) contributes with 13% and the order 0(e 2 ) with 6% to the value of 7 
(1.24), respectively. For the uniaxial Lifshitz case, the contributions for the 7^2 index are 
0.25 (17%) and 0.196 (14%). The very close values of the contributions of first and second 
order to 7^2 (as the €l parameter is 1.5 not being a small number), indicates that neglecting 
0(e|) could be a dangerous step in obtaining the exponent 7^2 via scaling relations in a 
more accurate way. Indeed, had we replaced the numerical values obtained for z/^ 2 = 0.73, 
r/L2 = 0.04 and d = 3 directly into the scaling law, we would have obtained 7L2 = 1.43. As 



argued in ]23| for the other critical exponents cyl2 and (3 12, whenever > 1 one should use 
the numerical values of z/^2, VL2 obtained from the e^-expansion for fixed values of (N, d, m) 
in order to obtain the numerical values of the other exponents via scaling laws. Therefore, 
provided we give this new interpretation to the numerical output of the e^-expansion when 
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tL > 1, we consider that the agreement between the numerical (Monte Carlo) and analytical 
(e^-expansion) results is remarkable for d = 3. The numerical value obtained here for the 
correlation length exponent is = 0.73. The experimental value of this critical index is 
still lacking. We hope our result sheds some light towards its experimental determination. 

The extension of the present method to the calculation of critical exponents for the 
m-fold (m 7^ 8) case reduces to the Ising-like case when m = and to the present case 
when m = 1 [pOR . An interesting open question is the calculation of the critical exponents 
Via and 77^4 using the e^-expansion at two- loop level. The approach followed here is not 
suitable to computing these critical exponents (parallel to the competition axis), since our 
choice of the symmetry point prevents a proper treatment in this direction. The possibility 
of devising another symmetry point to deal with these exponents seems to be feasible, and 
will be reported elsewhere. 



In recent articles, some authors [24] studied an alternative field-theoretic approach based 



on coordinate space calculations. In the first paper they recovered the results of reference [19 
for the cases m = 2, 6 analytically, but only could get the exponents numerically for the m = 
1 case, working entirely in coordinate space. It is worth emphasizing that they computed 
the fixed point only at one-loop order (see equation (82) in the mentioned paper). In the 
second paper, they computed the critical exponents at second order in perturbation theory 
by making use of a hybrid mechanism, going to coordinate or momentum space according 
to the necessity through a scaling function related to the free propagator in coordinate 
space. They obtained the exponents, whose coefficients of each power of are integrals to 
be performed numerically. The very similar values obtained for the exponents using their 
method or ours confirms that momentum and coordinate space calculations should give 
the same results, since either our approximation or the numerical approximation made by 
them is responsible for a rather small deviation in the two results when compared to 
the above numerical values. 

In conclusion, we have found a way to perform two- and three-loop integrals for the 
uniaxial Lifshitz point, needed to calculate universal properties at directions perpendicular 
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to the competition axis. The constraint in the loop momenta at the competition direction is 
the key ingredient to carry out the calculations. In this approximation, the loop momenta 
along the competition axis are not conserved when one uses this constraint. However the 
momenta along the (d — 1) directions are conserved. Momentum non-conservation along the 
competing direction as a higher-order effect does not seem to affect the critical exponents 
considered here in a significant way, as indicated by the comparison of our study with the 
available numerical data for the d = 3 case. 

It might be interesting to study general field theories including higher order derivative 
terms in this new framework. Topics including the extension of the present method to the 
region out of the Lifshitz point (S ^ 0) and two-loop calculations using a modified symmetry 
point along the competition axis are in development. 
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Figure captions 



• Figure 1. Feynman graphs corresponding to the integrals: (a) I 2 ; 
(b) I 3 ; (c) J 4 ; (d) I 5 . The broken lines in the graphs (b), (c), and (d) define the 
"internal" bubbles in each case. The momenta hi refer to the loop momenta in the 
quadratic and quartic directions, respectively. The labels pi, k[ denote the external 
momenta in the quadratic and quartic directions, respectively. 
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